Correlations of observables in chaotic states of macroscopic quantum systems 



o , 
o ■ 

(N ; 
^ ! 

Ph. 



> ■ 

(N 

O 

m 

O ' 

Oh: 



Ayumu Sug 

Department of Applied Physics, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka, 558-8585, Japan 

Akira Shimizi:0 

Department of Basic Science, University of Tokyo, 
3-8-1 Komaba, Tokyo 153-8902, Japan and PRESTO, 
Japan Science and Technology Corparation, 4-1-8 Honcho, Kawaguchi, Saitama, Japan 

(Dated: February 1, 2008) 

We study correlations of observables in energy eigenstates of chaotic systems of a large size A''. 
We show that the bipartite entanglement of two subsystems is quite strong, whereas macroscopic 
entanglement of the total system is absent. It is also found that correlations, either quantum or 
classical, among less than N/2 points are quite small. These results imply that chaotic states are 
stable. Invariance of these properties under local operations is also shown. 
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It is generally believed that almost all macroscopic sys- 
tems have chaotic dynamics, because otherwise thermo- 
dynamics would not hold. However, properties of chaotic 
quantum states in macroscopic systems are poorly under- 
stood as cornpared with those in systems of small degrees 
of freedom Among such properties, correlations of 
observables at different points are of particular interest. 
For example, such correlations are directly related to en- 
tanglement which is the central subject of quantum 
information theory The correlations also define the 
'cluster property', which is one of the most fundamen- 
tal notions of field theory Q. Moreover, the correlations 
determine the magnitudes of fluctuations of 'additive op- 
erators' (see below), which are macroscopic observables. 
Furthermore, Shimizu and Miyadera P| (Hereafter re- 
ferred to as SM) showed that two-point correlations de- 
termine the stabilities against classical noises, perturba- 
tions from environments, and local measurements. In 
particular, they showed that states with 'macroscopic en- 
tanglement' are unstable. Since chaotic dynamics is gen- 
erally believed to promote entanglement production [^, 
and since classical chaos is characterized by extreme sen- 
sitivity to the initial condition, it might be tempting to 
conjecture that chaotic states would be unstable. How- 
ever, as we will show, such a naive expectation is wrong. 
In this work, we study these points for chaotic quantum 
states in macroscopic systems, i.e., systems with a large 
but finite degrees of freedom. 

Two-point correlations and fluctuations of additive op- 
erators: We consider an energy eigenstate j^*) of a quan- 
tum chaotic system which is composed of N (^ 1) sites, 
where the Hilbert space is the tensor product of local 
Hilbert spaces. Let {aa(l)} be a basis of local observ- 
ables at a site /. Assuming that aaHYs are bounded, 
we normalize them as Tr[a]j(/)a^(/)] = constant x Sap- 



Then, all information about the two-point correlations 
in l^") are included in the variance-covariance matrix 
(VCM), Vai,0i' = {^\Saiil)Sap{l'm, where M„(/) = 
Oail) — (^'laa(OI^) This matrix also gives us in- 
formation about fluctuations (5'|(5ylt(5^|\E') of 'additive 
operators' A, which are defined as the sums of local op- 
erators; A = J2a i^aio-ail)- Here, Cq,/'s are c-numbers 
independent of N. Without loss of generality, we here 
normalize them as ^ \cai\^ = N. As shown in Ref. 0], 
the maximum fluctuation of the additive operators is 
Ne-axax, where emax is the maximum eigenvalue of the 
VCM. For example, Cmax = Oil) |0| for any product 
state I^E") — 0; IV'/), whereas Cmax = o{N) for 'vacuum 
states' of many-body physics in accordance with thermo- 
dynamics 0, 13] . Interestingly, as special states in large 
but flnite systems, there also exist pure states for which 
Gmax = 0{N) Q. Such purc states are superposi- 
tions of macroscopically distinct states, hence are called 
macroscopically entangled states 

Let us evaluate the VCMs of energy eigenstates in 
chaotic systems using random matrix theory (RMT). 
Suppose that an energy eigenstate |^'a), labeled by a 
single parameter A, is represented as \^x) = '^iCxiV) 
in some basis The ensemble averages, denoted by 

overline, of products of the coefficients can be calculated 
easily using RMT [lOj. For example, in the d x d GUE 
or GOE, 
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where q = and 1 for GUE and GOE, respectively. Up 
to the 4-th order, the other combinations which are not 
written above vanish identically. By using these formu- 
lae, we can evaluate the averages and variances of the 
matrix elements of VCMs. Since the right-hand sides of 



2 



Eqs. Q-® are independent of A, we will hereafter drop 
the label A. Furthermore, we will denote (4'a| • I^'a) sim- 
ply by (•). 

We here present results for the case where each site 
is a qubit, which is equivalent to a spin-i system. In 
this case, we can use the Pauli matrices {ux, cTy, az} as 
a basis of local observables Therefore, the VCMs 

of an A^-qubit system are 37V x 3A^ matrices. It is easy 
to show that {cra{l)) = and {aa{l)){(J(3{l')) — for 
I 7^ /' (a, f3 = x,y, z). With the help of these, the average 
values of the elements of the VCM are calculated as 
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This suggests that the VCM converges to the unit ma- 
trix / as — > oo, hence Cmax = 0{1). To examine this 
point, we plot in Fig. ^ the maximum and the minimum 
eigenvalues Cmax and eniin7 of the VCMs of energy eigen- 
states of CUE and GOE as functions of A^. It is found 
that both Bmax and emin (hence all eigenvalues) concen- 
trate on 1 as increases, and thus the VCMs approach 
/ .This means that fluctuations of all additive operators 
are 0{N) or less, and that all two-point correlations for 
I 7^ /' are vanishingly small. In fact, for I I' we can 
directly show that 
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hence \{5aa{l)Saf3{l'))\'^ is also 0(2^^). Therefore, either 
entanglement or classical correlation will be found to be 
exponentially small if one observes arbitrary two points 
of a chaotic state. Thus, in the sense of SM, chaotic 
states have the 'cluster property' [l5|, and are not entan- 
gled macroscopically. According to the general results 
of SM, these imply that chaotic states are stable under 
local measurements, and that the decoherence rate F of 
chaotic states is 0{N) or less under weak perturbations 
from any noises or environments which have short time 
correlations. This is in sharp contrast to macroscopically 
entangled states for which some noise gives F ~ 0{N'^) 
0. These points will be discussed later again. 

Correlations among many points: We have seen that 
two-point correlations are infinitesimal. How about m- 
point correlations for larger ml We can estimate it from 
bipartite entanglement between two subsystems as fol- 
lows. 

Fet us separate the system into two subsystems A and 
B which contain A'a and A^b = N — Na sites, respec- 
tively. The Hilbert space of subsystem A (B) has dimen- 
sion cJa = (de = 2^^). We assume that A^a < A^b, 
and evaluate 'purity' Tr {p\) as a measure of bipartite 
entanglement, where pA is the reduced density operator 
of A; PA = TrB [3 The purity takes the max- 

imum value 1 when the state |^) is separable (unentan- 
gled), and the minimum value I/cIa when pA is a scalar 
matrix I/cIa- For QUE and GOE, Eqs. O-© yield the 




FIG. 1: The maximum and the minimum eigenvalues of the 
VCMs in spin- 1 systems as functions of the system size A'^. 
The solid lines represent the results for eigenstates of GUE. 
The results from GOE, which are not displayed here, are very 
similar to those of GUE. The dotted lines are the results for 
the the central (2^~^-th) eigenstates of the hamiltonian ||HJ. 
In both cases the averages have been taken over 100 samples. 
The error bars show the standard deviation. 
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where AA^ = A'b — Na (> 0), and '• • •' denotes smaller 
terms. It is found that the average purity is quite small, 
which means that the bipartite entanglement of two sub- 
systems are strong. With increasing AA'^, in particular, 
the average purity approaches exponentially the mini- 
mum value 1/dA. Therefore pA converges to I/dA expo- 
nentially with increasing AA^, for almost all states. This 
is demonstrated for A'^ = 12 in Fig. 12 in which pm = Pa 
and m = A'a- It is noteworthy that the bipartite en- 
tanglement is large but not maximum when A'a = A^b- 
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Suppose 

diih), 0,2(12), 

hj2i'''Jm are measured. Since these observables 
commute with each other, their correlations can be 
calculated in a manner similar to the case of classical 
stochastic variables. In particular, the cumulants are 
given by 



{di{li) . . . drnikn)) c = (-1)'' 



d"' In Z 
dJi . . . dJ„ 



(7) 



where Z{{Ji\) = {^\ex:l>[-Y,^J^d^{li)\\'^) is the gen- 
erator of moments. If we regard the set of to sites 
^1,^2,' ■■j^m as subsystem A, we have Z{{Ji\) — 
Tr(p„exp[-^™^ Jiaj(?i)]), where pm = Pa- When 



N-m 



I/T 



3 





— Maximum X 


- 


GUE / 




GOE A::'^-''' 




Local int /^'^^ 

i—\J\J\A.\ II ILi /y ^ 


- 









1 2 3 4 5 6 

m 



FIG. 2: — logj Tr (pm) of a subsystem composed of m sites 
is plotted for N=12. The solid line denotes the upper limit, 
corresponding to the lower limit of Tr (pm), for each m. The 
dotted and dashed Unes, which look overlapped each other 
in this figure, represent results for GUE and GOE, respec- 
tively. The long dashed line shows the result obtained from 
the central (2^~^-th) eigenvector of the hamiltonian ((SJ. 



to a good approximation, according to our results 
above. In this case, we obtain {ai{li) . . .dm{lm))c — 

0, which means that there is no correlation among 
ai(/i), . . . , am{lm)- Since they are arbitrary observables 
at m points, it is concluded that neither entanglement 
nor classical correlation can be found if 4"^-^/^ < l, 

1. e., if the number m of points one performs measure- 
ments is somewhat smaller than half the total number of 
points. 

On the other hand, if observables at all points are mea- 
sured, one can detect all types of entanglement. For ex- 
ample, one can detect the bipartite entanglement by mea- 
suring a correlation of the Clauser-Horne-Shimony-Holt 
(CHSH) type This correlation becomes almost max- 
imum for chaotic states, because we have already shown 
that the bipartite entanglement between subsystems A 
and B is almost maximal (when 2'^^ ^1). However, the 
enhancement of the CHSH correlation over that of sepa- 
rable states (or of local classical theories) is at most by 
a factor of ^/2. Recently, another correlation of local ob- 
servables was found in Ref. , which detects macroscopic 
entanglement in general states (either pure or mixed): It 
becomes as large as 0{N^) for macroscopically entan- 
gled states, whereas it is at most 0{N) for states with 
Gmax = 0(1) (and for any mixtures of such states). For 
chaotic states, this correlation is at most 0{N) because 
we have shown that emax = 0{1). 

A model with short-range interactions: We have as- 
sumed RMT in the above analysis. One might suspect 
that RMT would not be fully applicable to macroscopic 
systems, in which interactions are short-ranged, because 



RMT assumes that strengths of interactions are of the 
same order between any two points. To explore this 
point, we study eigenstates of the following hamiltonian: 

N 

H = J^{a,(Z)a,(/ + l)+a,(Oa,(Z-|-l) 
(=1 

-f\/2cos 4)1 ay{l)ay{l + 1)} 

N 

-/i^{sin6l,a:r(0 +cos6';ct40}- (8) 

1=1 

Here, J and h are constants, {<j)i} and {9i\ are random 
variables with < (/i;,^; < 2tt, and cr^iN + 1) = cra(l). 
This hamiltonian describes a one-dimensional spin sys- 
tem in which spins interact only with their nearest neigh- 
bors, where the interaction in the y-direction is random. 
Besides, there is an external magnetic field with strength 
h, whose direction is random in the x-z plane. When J 
and h are sufficiently large (say, J — h — 1), this system 
becomes chaotic, except for states around the lower and 
upper limits of the energy spectrum, in the sense that the 
level spacing distribution coincides with that of GOE. 

The dotted lines in Fig. ^ represent the maximum and 
the minimum eig envalues of the VCM for the 2^-i-th 
eigenstate, which is located at the center of the spec- 
trum, where the chaotic nature appears most clearly. It 
is seen that the results agree fairly with those of RMT. 
Furth ermore, the long dashed lines in Fig. 13 represent 
— log2 Tr (/5^) for the same state. The results are close 
to those for RMT, especially when m is small. Simi- 
lar results are obtained also for other eigenstates except 
those near the upper and the lower ends of the spectrum. 
Since the density of states has a dominant peak at the 
center of the spectrum, our conclusions from RMT hold 
for the great majority of the eigenstates of the Hamilto- 
nian ((SIi.Tt'I We thus see that RMT correctly describes 
correlations in energy eigenstates of the chaotic system 
with short-range interactions. 

Invariance under local operations: What happens if 
local operations are performed on a chaotic state? Let 
{\ii)} be a basis of the local Hilbert space at site I. An 
energy eigenstate Ivp^^)) of an A^-site system can be ex- 
panded as 

i*'"^> = EE---E4"Uki)N2)---m- (9) 

ii 42 ijv 

For GUE, the probability density of the 2^ coefficients 

e^*(i:-i:ir„r-i)- (10) 

Since this density is invariant under changes of the lo- 
cal basis {\ii)} for any the statistical properties of 
|v[/(^)^ are not changed by local unitary transforma- 
tions. Furthermore, when a local projective measure- 
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ment that diagonalizes {|iAr)} is performed on the A'^- 
th qubit, the post-measurement state (for each measure- 
ment) is given by |?Ar) ■ ■ ■ Ki> ' ' ' Na") = 
1*^)1^'^^'^), where N' = N - 1, and c^^'l^, = 
cif.'.,„,.„/(E.,^-;-E.„, It is ea"sy to 
show that {cl^..\^,} also obeys the probabihty distri- 
bution of GUE, i.e., Eq. ((Tn|l with N being replaced 
with N'. Therefore, all the resuhs for |*(^)) hold for 
|*(^')) as weh, with N being replaced with TV' [Tl |. 
The same can be said when a local projective measure- 
ment that diagonalizes another local basis {Ijat)} is per- 
formed on the iV-th qubit, because the density is 
invariant under changes of the local basis. The same 
conclusions can be derived for GOE as well. We there- 
fore conclude that properties of chaotic states which we 
have found in this work are invariant under 'local oper- 
ations, ' including local unitary transformations and lo- 
cal projective measurements diagonalizing a local basis. 
This implies, for example, that a chaotic state cannot 
be disentangled by local measurements at less than N 
points. This is in sharp contrast to, e.g., a 'cat state,' 
|C) = I tt • • • T)/V2 + Ui • • • i)/V2, for which a single 
local measurement suffices to disentangle the state [3. 
This is consistent with the conclusion, which we have de- 
rived above using the theorem by SM, that chaotic states 
are stable under local measurements because two-point 
correlations are infinitesimal. 

Discussions: If one defines the 'degree' (or 'strength') 
£ of entanglement by the minimum number of local oper- 
ations that are necessary to disentangle it, then, accord- 
ing to our results, £ of chaotic states in large systems is 
quite high. However, this does not mean that they are 
anomalous as macroscopic states. Indeed, we have shown 
that they are stable against local measurements and that 
fluctuations of all additive operators are 0{N) or less, in 
accordance with thermodynamics Furthermore, 
as we have mentioned above, the correlation of Ref. 0, 



which detects macroscopic entanglement, is quite small 
for chaotic states. Moreover, in the infinite-size limit, 
all multi-point correlations vanish since the number m 
of measured points is always finite (hence m <C N/2 as 
N 00). Although the absence of quantum correlations 
for finite m at — > 00 has been generally proved |^, 
we have found here that chaotic states do not even have 
classical correlations in the same limit. Therefore, corre- 
lations, either quantum or classical, of chaotic states are 
neither detectable nor usable in infinite systems. 

N -dependences: In this work, we have often discussed 
N dependences of various quantities. Unlike a uniform 
state in a uniform system, however, correspondence be- 
tween states in systems of different sizes is non-trivial in 
chaotic systems. For completeness, we finally describe 
what the N dependences mean in this paper. 

For each system size N, consider a set S of all en- 
ergy eigenstates whose energies are in an interval {E — 
A/2, E A/2]. To see the N dependence of some quan- 
tity Q, look at its (probability) distribution P/v((5) in S. 
Then, the N dependence of Q can be discussed in terms 
of the N dependence of Pn{Q)- For example, we say 
'Q < Qo for almost all states in S for sufficiently large 
if for arbitrary positive number e there exists Nq such 
that 1 - Pn{Q < Qo) < e for N > Nq. We can thus 
discuss N dependences of various quantities for chaotic 
states. In particular, we can apply the results of SM, 
in which N dependences play crucial roles. Note that 
Pn{Q) is independent of E and A in the case of RMT. 
In the short-range interaction model given by Eq. ||SJ), 
the dependence of Pn{Q) on E and A is quite weak if 
we confine ourselves to states around the center of the 
energy spectrum. In either case, we can discuss the N 
dependence without specifying E and A. 
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